Pulse retrieval and soliton formation in a non-standard scheme for dynamic 
electromagnetically induced transparency 
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We examine in detail an alternative method of retrieving the information written into an atomic 
ensemble of three-level atoms using electromagnetically induced transparency. We find that the 
behavior of the retrieved pulse is strongly influenced by the relative collective atom-light coupling 
strengths of the two relevant transitions. When the collective atom-light coupling strength for the 
retrieval beam is the stronger of the two transitions, regeneration of the stored pulse is possible. 
Otherwise, we show the retrieval process can lead to creation of soliton-like pulses. 

PACS numbers: 42.50.Gy, 03.67.-a, 42.65.Tg 



I. INTRODUCTION 

Recent progress in the coherent control of light-matter 
interactions has led to many interesting possibilities and 
practical applications. Amongst them is the concept 
of electromagnetically induced transparency (EIT), first 
proposed by Harris jl|, in which a strong coherent field 
( "control" ) is used to make an otherwise opaque medium 
transparent near atomic resonance for a second weak 
("probe") field. This EIT scheme can be used to store 
and retrieve the full quantum information in a weak 
probe field by changing the strength of the control field 
while the pulse is inside the atomic sample [2(. In this 
paper we examine in detail an alternative method of re- 
trieving the stored information that was first proposed by 
Matsko et al. and investigate the parameter regime 
under which this process is feasible. 

The usual way of retrieving the stored information con- 
sists of a time-reversed version of the writing process 
0, . The elegant physics behind this scheme was de- 
scribed by Fleischhauer and Lukin, who noted that the 
combined atomic and optical state adiabatically follows 
a dark state polariton [6| . The "writing" process involves 
turning the control field to zero, storing the quantum in- 
formation of the light beam as the purely atomic form 
of the polariton. When the control field is returned back 
to its original value, the polariton switches back to pho- 
tonic form, identical to original optical pulse. In work 
describing some of the detailed behavior of that process, 
Matsko et al. noted an alternate scheme that may also 
store and reproduce a copy of a weak probe pulse 0. 
In this scheme, the writing process remains the same, 
but retrieval of the pulse involves applying a coherent 
control field to the transition originally coupled by the 
probe field. This causes the probe pulse to be regener- 
ated on the transition originally coupled by the control 
field. This alternative scheme cannot be explained in 
terms of dark state polaritons, and behaves quite differ- 
ently for different parameters of the system. We analyze 
this scheme in detail in this paper. 

As it is important to distinguish between the different 



control fields applied at different times, we will refer to 
the control field as the "writing" field during the first 
part of the process (storage) and as the "retrieval" field 
during the second part (the regeneration of the probe 
pulse). 

We will show that when the collective atom-light cou- 
pling strength for the retrieval beam is the largest of 
the two transitions, retrieval is possible and the retrieved 
pulse is amplified, time reversed, stretched or compressed 
in time and phase conjugated compared to the input 
pulse. Conversely, if the collective coupling strength of 
the retrieval beam is not the largest of the two transi- 
tions, we find that the retrieved pulse differs substantially 
from the input pulse and the retrieval process can lead 
to the creation of a soliton-like combination of electro- 
magnetic fields and atomic coherences that propagates 
without change in shape. 



II. MODEL 

To analyze the system we use a quasi one-dimensional 
model, consisting of two copropagating pulses passing 
through an optically thick medium of length I consist- 
ing of three-level atoms. The atoms have two metastablc 
ground states \b) and |c) which interact with two fields 
£ p (z,t) and £ c (z,t) as shown in Figure ^ &i are the 
slowly varying amplitude related to the positive fre- 
quency part of the electric field given by 



Here lo^ = [E^ — E v )/T% is the resonant frequency of 
the \fi) <-> \v) transition and V the quantization volume, 
here taken as the interaction volume. As the pulses are 
co-propagating we are able to neglect Doppler effects. 

To perform a quantum analysis of the light-matter in- 
teraction it is useful to use locally-averaged atomic oper- 
ators. We take a length interval Sz over which the slowly- 
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a small input pulse of maximum amplitude f2° <C £1° is 
sent into the medium on the \b) <-» \a) transition. EIT 
effects cause the group velocity of the the input_pulse to 
be drastically reduced to a new value given by 
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FIG. 1: Level structure of the atoms 



varying field amplitudes do not change much, containing 
nASz 3> I atoms, where n is the atomic density and A is 
the cross sectional area of the pulses, and introduce the 
continuous atomic operators 
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where Zj is the position of the ith atom and \p l (t)} is the 
state wavefunction for the ith atom.. 
Using these continuous atomic operators, the interac- 
tion Hamiltonian under the rotating wave approximation 
is given by 
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Nh 



[g p £ p {z, t)a ab (z, t)+g c £ c (z, t)a ac (z, i)+H.c. 

(2) 

where I is the length of the cell, N is the number of 
atoms in the interaction region and the coupling con- 
stants are g p = d ab ^ujab/2e Vh, g c = d ac ^LO a c/2e Vh 
where d a b and d ac are the dipole moments of the \a) <-> \b) 
and | a) <-> \c) transitions respectively. This Hamiltonian 
leads to the following equations of motion for the density 
matrix elements and fields 
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which can be much less than the speed of light. This leads 
to significant spatial compression of the input pulse as it 
enters the medium, allowing the entire pulse (the typical 
pulse length outside the cell is of order of a few kilome- 
ters) to be stored inside a vapor cell (of length a few 
centimeters). Once the input pulse is completely inside 
the medium, we slowly turn off the writing beam on the 
|c) «-» | a) transition. This writes the information carried 
by the input pulse onto a collective atomic coherence for 
storage @. 

After a controllable storage time T s , the stored infor- 
mation can be retrieved by turning on a retrieval beam 
driving the \b) <-> \a) transition as proposed in 0, Q. 
Note that this is in contrast to the usual retrieval scheme 
for EIT light storage in which the retrieval beam is on the 
\c) <-> | a) transition and the output pulse is generated on 
the 1 6} «-» | a) transition |(J. For normal EIT in the ideal 
case, the output pulse is identical to the input pulse. 
Here, because the time reversed version of the writing 
I-beam [ s on ^ ne |^ <_> transition, the output pulse is 
generated on the |c) <-> \a) transition. As a consequence 
the output pulse need not have the same frequency or 
polarization as the input pulse. 

It is not immediately clear why any fields made by 
this process would be correlated with the stored pulse. 
After all, the original pulse produces a small population 
in state \c) and leaves most atoms in state |6), so the 
control beam is mainly interacting with states that have 
not been affected by the original pulse in any significant 
way. Also, as the strong retrieval beam is initiated on 
the most populated transition, it will lead to significant 
spontaneous emission, which would appear to dominate 
any kind of coherent process necessary for the retrieval. 
However, as will be shown later, our simulation indicates 
that under some parameter regimes, it is indeed possible 
to recover an output related in amplitude and phase to 
the input pulse. 



III. RETRIEVING THE PROBE PULSE 



where the 7 M are phcnomcnological decay rates and we 
have defined the two Rabi frequencies as f2j = gi{£i) and 
the constants a p = N\g p \ 2 and a c = N\g c \ 2 are the col- 
lective atom-light coupling constants for the transitions 
1 6) <->■ | a) and \c) «-» \a) respectively. 

We now proceed to consider the modified version 0, 
of the storage and retrieval scheme using dynamic EIT. 
With all the atoms initially optically pumped into the 
ground state |6), we turn on the writing beam driving 
the |c) <-> | a) transition to its maximum value 17", while 



We will start with an example that illustrates a suc- 
cessful pulse retrieval, and then go on to analyze other 
possible behaviors of this system. We assume standard 
atomic initial conditions for EIT given by pbb(z,0) = 1; 
p cc (z,0) = pab(z,0) = p c b(z,0) = p ca {z,0) = 0. In or- 
der to better identify the properties of the new retrieval 
scheme and to judge the quality of the retrieval process, 
we choose an input pulse at z = to be the sum of two 
Gaussians of different height (so the total input pulse is 
non-symmetric) with a time-varying phase factor. The 
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FIG. 2: Envelope of pulses entering the medium at z — as 
a function of time, a) Q p is the field driving the \b) <-» \a) 
transition and b) £l c is the field driving the |c) <-> \a) tran- 
sition. The parameters are T q jj = 1407 _1 , T s = 18.857 -1 , 

T on = 2597 -1 , Q° p = 0.02657 and = 2.65267- Amplitude 
of pulses are displayed in units of 7. The input pulse to be 
stored is Q p before the start of storage. The retrieved pulse 
is Q c after the storage process. 



FIG. 3: Comparison of a) amplitude and b) phase of both 
the retrieved and input pulses. The solid line represents the 
retrieved pulse and the dashed line the input pulse (magnified 
by a factor of ten). The parameters are a p = 301777, a c = 
29272C7, I = 4cm. The input pulses are as shown in Figure 
121 The large amplitude on the retrieval pulse transition at 
earlier times is because the field on this transition was used 
as the writing beam initially. Similarly the large value on the 
input pulse transition at later times represents the retrieval 
beam. 



boundary conditions for the fields at z = is of the form 



n p (o,t) = n£/(t)e*w + 
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where f(t) is a unit amplitude envelope function. The 
first part of equation JHJ) describes the input (signal) pulse 
that we wish to store, whereas the second part of equa- 
tion © describes the turning on of the retrieval beam 
(amplitude Q®) at time T on with a switching time of ap- 
proximately T s . fl c describes the turning off of the writ- 
ing beam at time T q jj. These boundary conditions for 

the fields are plotted in Figure [21 

With these initial conditions, the equations of motion 
ijljl can be solved numerically in the moving frame defined 
by£ = z, t = t — z/c using the method described in 
Shore implementing the numerical integration with a 
fourth order Runge-Kutta method 0. For convenience, 
we choose 7t = -y a b = 7 an d 7 C = 7ca = (a c /a p )j. 

Figure |21 shows a comparison of the phase and ampli- 
tude of the retrieved and input pulses, both propagat- 
ing in the positive z direction. In practice, this output 
pulse of well-defined shape will be superimposed on top 
of spontaneously emitted photons, so for the purposes of 
observing this output, it would be better to choose the 



lowest atomic density allowed for EIT to work. 

From the graphs, we observe that the retrieved pulse 
is time reversed, amplified, widened in time and is the 
phase conjugate of the input pulse. We can understand 
this behavior by examining the interaction of the retrieval 
beam with the stored coherence and how the retrieved 
pulse is generated. As the writing part of our process is 
identical to the usual dynamic EIT setup @, we know 
that at the end of the storage process, the only variable 
of the system aside from p^j, that is significantly nonzero 
is p c b whose spatial variation, shown in Figure^ encodes 
the phase and amplitude information of the original input 
pulse. 

The retrieval beam enters the medium driving the tran- 
sition I b) <-> I a) while almost all the atoms are in state |6). 
As the medium is optically thick, the retrieval beam is 
strongly absorbed and its initial wavefront moves across 
the medium at a speed much less than the speed of light. 
The spatial variation of the retrieval beam at any time 
during the retrieval process has the general shape shown 
in Figure [SI At any time, this shape divides the medium 
into three regions with distinct dynamics. In region I, the 
retrieval beam Q p has attained its maximum value as it 
has optically pumped all atoms into the state |c). In re- 
gion II, the wavefront of the retrieval beam is severely 
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FIG. 4: Amplitude a) and phase b) of p c b during storage. 
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FIG. 5: Spatial variation of the retrieval beam divides the 
medium into three regions as described in the text. At any 
time, the retrieved pulse is generated due to the dynamics 
that occur in region II where the retrieval beam is in the 
process of pumping population out of state \b). In region I, 
all atoms are in state c). The beam has not yet penetrated 
to region III. 



attenuated due to absorption and in region III, not yet 
reached by the retrieval beam, all the atoms are still 
in the polariton state left by the writing process. This 
means that most of the atoms are in state \b) here. We 
note that at any point in time, the interesting dynam- 
ics related to the retrieval of the stored information is 
occurring only in region II. In this region the retrieval 
beam f2 p is in the process of pumping atoms from \b) to 
a) at a rate of \fl p \ 2 /^. Before a significant population 
has accumulated in |c), the retrieval beam can coherently 
scatter off the stored coherence p c b, generating a retrieved 
pulse on the \a) to \c) transition. These generated pho- 
tons will encode some of the properties of the input pulse 
and propagate out of the medium at the speed of light 
through region III, as this region contains only atoms in 
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FIG. 6: a) p cc (dashed) and p aa (solid) as a function of 
z during the retrieval process. A significant population re- 
sides in the excited state during the retrieval process, b) 
|p c &(0.024, t)\ as a function of tj" 1 . The spike in the graph 
shows the effect of optical pumping by the retrieval beam 
which overwrites the stored coherence. 



On long time-scales optical pumping alters p c b and will 
overwrite any coherence that was initially stored there. 
At this point, photons emitted when atoms decay from 
\a) to |c) will bear no relation to the input pulse. Further- 
more, any field on the |c) <-> |o) transition will be strongly 
damped due to the significant population in state \c). 

Note that in contrast to the adiabatic following in the 
usual EIT procedure .6] , the retrieval process here is non- 
adiabatic and as demonstrated in Figure EJa), the pop- 
ulation in the excited state |o) can be quite significant 
during the retrieval process. 

Using this description of the dynamics, one can explain 
many of the features exhibited by the retrieved pulse. 
The time reversal occurs because the wavefront of the 
retrieval beam moves slowly from left to right, so the part 
of the stored coherence that is closest to the cell entrance 
sees the retrieval beam first. This part of the coherence 
near the medium entrance corresponds to the tail end 
of the input pulse (i.e. the part of the input pulse that 
entered the cell last), so the tail end of the input pulse 
will be retrieved first, resulting in an output that is time 
reversed. Figure |3| also demonstrates that the retrieved 
pulse can be amplified, and in this particular example 
the amplitude is increased by a factor of about twenty. 
This occurs because the new retrieval scheme has a larger 
reservoir of atoms available for producing photons. In the 
normal EIT setup photon number in the retrieved pulse 
is limited by the number of atoms in state |c) during the 
storage process which is strictly less than the number of 
photons in the input pulse. In this scheme the number of 
photons generated is limited by the maximum number of 



5 




50 100 150 200 250 300 350 




1 50 100 150 200 250 300 350 

FIG. 7: Comparison of a) amplitude and b) phase of retrieved 
(solid) and input (dashed) pulses. The input pulse has been 
magnified by a factor of ten. The amplitude of the retrieval 
beam is now = 5.30527 twice the value in Figure[3] We see 
that the retrieved pulse is made earlier in time and narrower 
in time because of the increased pumping rate by the retrieval 
beam. All other parameters are identical to those in Figure 

m 




FIG. 8: Comparison of a) amplitude and b) phase of retrieved 
(solid) and input (dashed) pulses. The input pulse has been 
magnified by a factor of ten. The abrupt cut off of the re- 
trieved pulse near rj" 1 — 600 is due to the arrival of the 
retrieval beam at the exit of the cell z = 0.04 which pumps 
all atoms into state |c) and prevents any field existing on this 
transition. The parameters are a p = a c — 301777, I = 4cm. 
The input pulses are as given in Figure |5] 



atoms in state c) that the medium can sustain without 
destroying the stored coherence. This is a property of 
both the size of the stored coherence and the medium 
(e.g. optical density). It is, however, independent of the 
properties of the retrieval beam. For an optically dense 
medium, this is generally much larger than the number 
of photons in the input pulse. 

It is also clear that the width of the retrieved pulse is 
determined by the speed at which the wavefront of the 
retrieval beam propagates across the medium and not the 
initial pulse length. Figure \7\ shows the retrieved pulse 
when the amplitude of the retrieval beam is doubled com- 
pared to Figure[3]while keeping all other parameters iden- 
tical. We see that the output pulse is narrower and its 
amplitude larger while the total number of photons (as 
indicated by the area underneath the graph of intensi- 
ties) remains approximately the same. This is due to the 
the more intense retrieval beam now being able to pump 
atoms out of \b) at a higher rate, enabling it to move 
across the medium more quickly and ensuring a shorter 
time interval between the retrieval of the front and back 
part of the original input pulse. Since the total number 
of generated photons is unaffected by the intensity of the 
retrieval beam, energy conservation necessitates that the 
retrieved pulse has greater amplitude. 

The phase conjugation of the output pulse is most eas- 
ily understood by examining the equations of motion Q 
and ©. On a short time scale, we have pbb ~ 1 S> 
PccPaa, and the following set of equations describes the 
generation of the retrieved pulse from the conjugate of 



the stored coherence p* cb 

[at + c q- z )^ = w ^ ( g ) 

Pac (z,t) = i [ e-^ t -^Cl p {z,s)p* cb {z,s)ds. (9) 
Jt ok 

Since p* b stores the conjugate of the input phase, the 
retrieved pulse is therefore the phase conjugate of the 
input pulse. 

From examining the phase of the generated output field 
compared to the input, it is clear that the quality of the 
retrieval process, (in terms of extracting a pulse of the 
same shape as the time reversed input) is not perfect. 
For example, from Figure (b), we see that the phase of 
the retrieved pulse fails to decay to zero after a certain 
time. This feature is even more prominent when we set 
the collective atom-light coupling constants to be equal 
a p = a c , as shown in Figure |H1 Again we see that the 
phase of the retrieved pulse plateaus, this time near the 
peak of the second Gaussian. Increasing a c further to 
move into the regime a c > a p , we see from Figure [5] 
that the quality of the retrieval process is so bad that 
the output pulse does not even display the characteristic 
double peak of the input pulse. 

To further investigate the behavior of the retrieved 
pulse as the relative strength of the collective coupling 
constant is varied, we consider a Gaussian input enve- 
lope as this makes clearer the difference in retrieval qual- 
ity. Figure ITUI show the three cases a p > a e , a p = a c , 
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FIG. 9: Comparison of a) amplitude and b) phase of re- 
trieved (solid) and input (dashed) pulses. The input pulse 
has been magnified by a factor of ten. The parameters are 
a p = 30177C7 j ct c — 31082C7, I — 4cm. The input pulses are 
as given in Figure |5] 
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FIG. 10: Comparison of amplitudes of input pulse (dashed) 
magnified by a factor of ten and output pulse (solid) for the 
three parameter regimes a p > a c (a), a p — a c (b) and a p < 



and a p < a c . For a p > a c , the generation of the out- 
put pulse automatically ceases so that the retrieved pulse 
will always display the falling end of the Gaussian input. 
However, as a c approaches the value of a p , the pulse 
generation takes a longer time to cease giving the assy- 
metric long tail shape shown in Figure EH (a) . When 
dp = eventually a constant value of the generated 
field is maintained until the retrieval beam reaches the 



exit of the cell and pumps all the atoms to |c). Continu- 
ing the trend, Figure ITUl (c) shows that in the regime of 
«p < a c , the generated pulse amplitude grows until the 
retrieval beam has optically pumped all the atomic popu- 
lation into |c). The phase variation of the retrieved pulse 
confirms the general trend that quality of the retrieval 
process decreases as we move across the three parameter 
regimes. 



IV. SOLITON-LIKE BEHAVIOR 

Our numerical results strongly indicate the existence 
of a soliton solution at the critical point a p — a c . At this 
point the field fl c generated from the stored coherence is 
able to induce a nonzero coherence p c b by a co-operative 
action with the retrieval beam Q p . The induced p c b in 
turn can then be used to generate £l c . This steady state 
cycling action explains why the system can continue to 
generate an output even when the 'dynamic' region (re- 
gion II in Figure |3J) has moved beyond where the in- 
put pulse was originally stored. More specifically, when 
«p > ct c , the induced coherence is not large enough to 
maintain the cycling action, causing the effect to die out; 
if a p = a c , the effect is self sustaining where the induced 
p c b is just sufficient to maintain the value of f2 c that gen- 
erated it; and for a p < a c , a greater coherence (p c b) 
is induced, leading to an output at the cell exit that is 
continually amplified in time until cut off by the arrival 
of the retrieval beam. However, we have observed that 
when the cell length is increased, the amplitude of the 
output pulse appears to tend to a limiting value. 

Figure ITU shows the induced p c b-, fl p and Q c for a p — 
a c when the initially stored coherence no longer exists. 
We note that the shape and size of p c b and the fields 
propagate unchanged across the cell, indicating behavior 
characteristic of solitons. We also found that the shape 
and size of the final output is independent of the values of 
the coherences originally stored, indicating that the final 
field-coherence formation is a characteristic of the system 
independent of the storage/retrieval process requiring an 
initial nonzero pb c coherence only as a seed. 

To demonstrate analytically that our results are in fact 
solitons we substitute the ansatz Qi(z,t) = Sli(z — vt), 
Ptiu(z, t) — p^(z-vt) into equations (QJ for the case a p = 
a c , where v is the soliton parameter that designates the 
constant speed with which the soliton propagates across 
the medium. This allows us to determine the following 
relationship between the two pulses and the coherence 
generated as 



Pcb(s) 



OtQV 



(10) 



where a p = a c = ao is the collective light-atom coupling 
constant identical for both transitions and s = z — vt. 
We also obtain the following relationship between the 
limiting values of the two fields 
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FIG. 11: |n p | (dash), |fi c | (dash-dot) and \p c b\ x 5 (solid) as 
a function of z at times t^ 1 = 477.6 (a), 489 (b), 501.6 (c) 
for a p = a c case. The field and atomic coherence propagates 
together unchanged across the cell 



a v 



(|fi2T + |fi?T) 



(ii) 



where f2^° = lim s ^_oo fij»(s) and f2£° = lim^oo Cl c (s). 
When we compare equations (|10f> and against the re- 
sults of our numerics for various soliton parameters (the 
soliton parameter can be changed by varying the am- 
plitude of the retrieval beam), we find good agreement. 
After further elimination of all the atomic variables, we 
obtain the two soliton equations relating the two fields 



n 



d 2 n* 

' ds 2 



ds 2 



9. 



dot 

ds 



c ds 



(12) 



d 3 n p 

'lb? 



d 2 n % 

ds 2 



dn 

ds 
2D 2 , d 
v ds 



d 2 n p 
'~dl 2 ~ 



dfl p 
i — - 

ds 



(\n c \ 2 + \n p \ 2 ) 



+ ^n}{\o?\*-\sic\ J -\n P \ J ) (is) 

from which the dispersive and nonlinear terms are clearly 
visible. 

We believe the photons generated from the cycling ac- 
tion outlined above bear little relation to the input pulse 
and are therefore not useful as far as quantum informa- 
tion retrieval is concerned. However the existence of the 
soliton-like solution and its sensitivity with respect to 
the coupling parameters is likely to lead to other inter- 
esting possibilities. A further analysis of the generation 
and properties of the solitons is beyond the scope of this 
paper, although the possibility that solitons can exist in 
atomic A systems has previously been considered 0, |nj 

In conclusion, we have demonstrated the feasibility of 
an alternative retrieval scheme for dynamic EIT under 
certain parameter regimes and provided physical expla- 
nations for its behavior. Our numerical simulation also 
demonstrated the ability of this new scheme to create 
soliton-like features that are sensitive to the relative cou- 
pling strength of the two transitions. Due to its sensitiv- 
ity to parameter change, we believe this solitonic behav- 
ior could prove useful within the context of magnetome- 
try or high precision measurement. 

We thank Elena Ostrovskaya for helpful discussions re- 
garding the soliton features. 
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